On the Decomposition of Lattice- 
Periodic Functions 
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(Manuscript received March 31, 1965) 

The problem of decomposing an arbitrary periodic function defined on 
an n-dimensional cubic lattice into finite linear combinations of certain 
primitive functions is considered. Generally, a primitive function is one 
which periodically assumes only the values ±1 and 0. Rather simple neces- 
sary and sufficient conditions are derived for such a decomposition and 
when a decomposition is possible, an algorithm is given which accom- 
plishes it. These results have been iised in recent generalizations of the Ewald 
method. 

I. INTRODUCTION 

In a study of the classical problem of the calculation of the potential 
due to an ionic crystal lattice, and in particular, generalizations of the 
Ewald method (Ref. 1) along the lines of Nijboer and Dewette (Refs. 
3, 4), W. J. C. Grant (Ref. 2) proposed the following problem: Suppose 
we say that an ionic crystal lattice is primitive if for a suitable choice of 
origin there exist three vectors .f i , .f 2 , f 3 such that the charge at the 
point rtiXi + n 2 .f 2 + Wa^a ' s i us ^ <?o ( — 1 ) " J n,+n, for some fixed g and for 
all triples of integers (wi , n 2 , n 3 ) and that the charge at all other points 
is zero. (For example, the ordinary NaCl lattice is primitive with the 
.t, taken to be the unit coordinate vectors and g = 1.) The question is 
then: Which crystal lattices can be decomposed into finite sums of 
primitive lattices? Different primitive lattices in the decomposition may 
have different origins and by the sum of two lattices we mean, of course, 
the component-wise sum. 

The object of this paper is threefold: 

(i) The problem is extended to its natural n-dimensional analogue. 

(ii) Rather simple necessary and sufficient conditions are given for 
the existence of the desired decomposition. 
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(m) When such a decomposition is possible, an algorithm is given 
which accomplishes it. 

II. PRELIMINARY IDEAS 

In order to illustrate the basic ideas which will be used in the proofs 
of the general (n-dimensional) theorem (seep. 1200), we begin by con- 
sidering the following one-dimensional version. 

Suppose we call a real-valued function / denned on the integers 
primitive if for some integers x and c it is true that 

./ v _ f(— 1)" if z =- ax + c 
— |p otherwise 

for all integers a. The question then becomes: What is the set of all 
those functions which can be represented as real finite linear combina- 
tions* of primitive functions? For example, the function g denned by: 



g(z) = 



if 2 = (mod 4) 

1 if 2=1 (mod 4) 

2 if 2 = 2 (mod 4) 
[-3 if 2 = 3 (mod 4) 



may be decomposed into a linear combination of primitive functions 
by: 



g(z) = -<7i (z) + 2^2 (z) + 03 (z) 



where 



f(-l)° if z = 2a 
9i[Z) "' \0 otherwise 

_ f(-l)° if z = 2a + 1 
9i{z) " \0 otherwise 

<7s (z) = (-I) 2 . . ■ 

We can write this more graphically if we usi the notation 

/: ■ • • , a , Oi , o 2 , • • ■ 

to denote the fact that /(0) = a ,/(l) = fli , etc. We then have 



* In this paper, linear combination will always mean^m7e linear combination. 
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-01 

9i 



, -1, 0,1, 0, 
, 0, 2,0,-2, 
, 1,-1,1-1, 



g: ... , 0, 1,2,-3, 

Similarly, if we start with 

*: ... ,1,3,-2,-1,-3,2, .-. 

then the desired decomposition is easily found to be: 

h: ..■ ,1,0, 0,-1, 0,0, ... 

3/i 2 : • • ,0,3, 0, 0,-3,0, ••• 

-2h 3 : ... ,0,0,-2, 0, 0,2, ••• 



fc x : 



,1,3,-2,-1,-3,2, 



In general, it is clear that any linear combination / of primitive func- 
tions is periodic and that within a period the sum of the function values of 
f must be zero. If a function has these latter two properties, we say that 
the function has mean zero. It might at first be surmised that any func- 
tion with mean zero could be written as a linear combination of primi- 
tive functions. However, a\ tempts to decompose the periodic function 



9- 



,1,-1,0,1,-1,0, 



(the bar indicating a complete period) soon lead one to suspect that 
this initial guess is incorrect. (In fact, g cannot be decomposed into 
primitive functions.) 

One question which arises immediately is exactly which periods the 
primitive components of a function / might have, if / itself has some 
period p (where we say that / has period p if f(z + p) = f(z) for all 
z). In the preceding example, while g has period 3, perhaps there is a 
decomposition of g for which the primitive components have much 
larger periods. (It will turn out, however, that this is not possible.) 

To answer these questions, we first introduce some notation. If g is 
a function defined on the integers,* then by g(z/r) we mean the func- 
tion defined by: 

g ( - ) if - is an integer 
otherwise 



* In general, in this paper all functions assume the value on points with non- 
integral coordinates. 
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Let i(z) denote the function which assumes the value 1 on all integers. 
Thus, if the function / which we wish to decompose has period 

p = 2 a (2m+ 1), 

then by forming the functions i[ (z — k)/ (2m + 1 )]/(z), ^ k < 2m + 1, 
we have functions which "sample" / at points separated by a distance 
of 2m + 1. For example, if/ is given by 



/: • • • , 1,3,-6,3,-5,4,1,3,-6,3,-5,4, • • • 

so that the period of / is 6 = 2-3 (where we will assume that /(0) = 1 ) 
then we have 



*(JJ/(z): •••,1,0, 0,3, 0,0, ... 
i (iZ-iW) : . . -,0,3, 0,0,-5,0, • • • 

i ( L ^)/0) : ' -,0,0,-6,0, 0,4, • • • 
Note that i[(z — k)/ {2m + 1 )]/(«) also has period p and, in general, 

The result toward which the remainder of this section will be devoted 
can now be expressed simply in the following way : If / has period 

p = 2 a (2m+ 1) 

then / can be expressed as a linear combination of primitive functions if 
and only if for each k the function i[(z — k)/ (2m + l)]f(z) has mean 
zero. 

It follows from this, for example, that if p = 2° then / can be decom- 
posed into primitive functions if/ has mean zero. On the other hand, if 
/ has an odd period p = 2m + 1, then each function i[ (z — k)/ (2m + 1 )] 
has just one nonzero value per period so that / can be decomposed into 
primitive function if it is identically zero. 

We now give a series of lemmas, informal proofs and examples which 
will indicate the ideas needed for the proof of the general theorem. An 
outline of our plan of attack is to establish the following results: 
If / is a linear combination of primitive functions then for 
any k and for any r ^ 0, f[(z — k)/r] also is a linear combi- 
nation of primitive functions. (1 ) 
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If / is a linear combination of primitive functions and / has 

period p = 2" (2m + 1 ) then for all k, i[ («'-*)/ (2m + 1 )]/(z) , * 

has mean zero. 

If / has period 2° and mean zero then / is a linear combination . . 

of primitive functions. 
Assuming we have established (1), (2) and (3), the proof of the 
original assertion follows directly. One direction follows immediately 
from (2). To show the other direction assume that for each k, i[(z — k)/ 
(2m -|- l)]/(z) has mean zero. Notice that each function i[(z — k)/ 
(2m + l)]/(z) is just an "expanded" copy of a function fk(z) which 
has period 2" and mean zero (i.e., i[(z — k)/(2m + l)]/(z) = fk[(z — k)/ 
(2m + 1)]). Hence, by (3), /*• (z) is a linear combination of primitive 
functions and it then follows by (1) that this is also true of /*[(z — k)/ 
(2m +1)]. Consequently 

f(z) =£i(^)/« = E/»(^) 

,i=o \2m + 1/ i-=o \2m + 1/ 

is a linear combination of primitive functions and the proof is completed. 

It remains to prove (1), (2) and (3). 

The proof of (1 ) is straightforward. We first note that if f(z) is primi- 
tive then f[(z — k)/r] is also primitive for any k and for any r ^ 0. 
For by hypothesis there exist x and c such that 

(—1)" if z = ax + c 



/U) " \0 otherwise 
On the other hand, by definition we have 

g _ A \f(y) if z = ry + k 



J 
Hence 



otherwise 



(^h: 



(-1)" if z = r(ax + c) + k 

otherwise 
(-1)° if z = a(rx) + (re + k) 
otherwise 



and so/[(z — k)/r] is primitive. The extension to linear combinations of 
primitive functions follows at once and (1 ) is proved. 

In order to prove (2) we first need an auxiliary result (a simplified 
version of Lemma 2 ) . This is : Suppose / has periods p = 2" (2m + 1 ) 
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and p = 2 a ' (2m' +1). Then for any k, i[(z - k)/(2m + 1 )]/(*) has 
mean zero iff i[(z — k)/(2m + l)]f(z) has mean zero. To prove this, 
let us first assume that a = a and 2m + 1 divides 2m + 1. The sum 
^f=o J i[(z — k)/(2m + 1 )]f(z) is exactly the sum of the f(z) for which 
z — k = (mod 2m + 1), and ^ z — k ^ p — 1 (since *[(« — fc)/ 
(2m + 1 )]/(*) has period p). There are 2" such values of z — k, namely, 

z - keA = {0,2m + l,2(2m + 1), ••• , (2 a - l)(2m + 1){. 

Similarly the sum 2*-° * *[(* ~~ W ( 2m ' + 1 )]/( 2; ) i s exactly the sum of 
the f(z) for which z — k = (mod 2m + 1 ) and ^ z— k ^p — 1. 
Again there are 2 a such values of z — k, namely, 

z- keB = {0,2m + 1, 2 (2m' + 1), ••• , (2 a - l)(2m' + 1)}. 

All the elements of A and B are congruent to zero modulo 2m + 1 
(since 2m + 1 divides 2m' + 1). Also since 2m + 1 and 2m' -f- 1 are 
odd then both sets A and B contain a complete residue system modulo 2 a . 
Hence modulo p, A and B are identical. Since / has period p then 

t=o \2m + 1/ zca zeB z=o \2m + 1/ 

If we now assume that a ss a (instead of a = a) then it is not diffi- 
cult to see that 

Thus, what we have shown is that if/ has periods p = 2 a (2m + 1) and 
p = 2"' (2m + 1) where p divides p then *'[(* ~ k)/(2m + 1)1/(2) has 
mean zero iff i[(z — k)/(2m + l)]/(z) has mean zero. Since in general 
a function which has periods q and q also has period (#,#' ) (the greatest 
common divisor of q and q), then the initial assertion follows at once. 
As a simple example consider the function / given by 



/: ••• ,1,3,-2,-1,4,2,1,3,-2,-1,4,2, 

This function has 6 = 2-3 as a period and i(z/S)f(z) has mean zero 
since 

£,i(z/Z)m = /(0) +/(3) = 1 - 1 =0. 

z=0 

However we may also consider / as having a period of 12 = 2 -3 in 
which case i(z/3)f(z) has also mean zero since 
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£t(«/3)/(«) = /(0) +/(3) +/(6) 4-/(9) = 0. 

2-0 

Finally,/ has a period of 18 = 2-3 2 and i(z/9)f(z) still has mean zero 
since 

X>-(z/9)/( 2 ) = /(0) 4-/(9) = 0. 

2-0 

Our next step will be to prove (2) using the result just established. 
We first show that if / is primitive and has period p = 2" (2m + 1 ) then 

E * f^MVVw = ° for a11 L (4) 

z =o \2ra + 1/ 

To see this, we partition the integers into two-element subsets {w t - , «,} 
such that Vi = u t , + 2m + 1 for each i. Since / is primitive there exist 
integers .t and c such that 

f(z) = ((- 1 )" if z = ax + c . 
\0 otherwise 

Since 

UiX = ViX (mod 2m -\- 1 ) 

and Vi — u, : = 2m + 1 is odd then it follows that 

f(uiX + c) = —f(viX -f c) for all i. 

But 

/gs - A = /g - fc\ and fc 

\2m + 1/ \2wi + 1/ 

Consequently it follows from the fact that / has period p that 

g'ten)'"- foralU: 

and (4) is established. 

To establish (2 ) assume that / has period p = 2 a (2m + 1 ) and is a 
linear combination of primitive functions /, , 1 ^ i ^ f. If /,■ has period 
p, = 2 n, (2m< 4- 1) then by (4) we know that 



S'feriV'W" forall,c - 
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Hence, if we choose q = pip 2 ■ • ■ p t p = 2 a (2m + 1) then by the sim- 
plified version of Lemma 2, we have 

<T i ( *~ k \fi(z) = for 1 £ i g t and all k. 

z =o \2m' + 1/ 

Consequently 

fc? / . _ Z- \ 

l/fe) = 



V • ( z ~ k 
h l \2m' + 1 



since by hypothesis / is a linear combination of the /, . But / has period 
p so applying the Lemma 2 result again we find 



S % \2m + l) 



f(z) = 



and (2) is proved. 

We are left with (3) to prove. To do this we first establish the follow- 
ing result: If h(z) is defined by h(z) = (—1)' then for a fixed n, the 
2 n - 1 functions h[(z - k)/2 T ], ^ k < 2 r , ^ r < n, are linearly 
independent over the reals. This is easy to see since for a fixed r, 
h[(z — k)/2 r ] assigns a nonzero value only to those z such that z = k 
(mod 2 T ). Hence for k = 0, 1, • • • , 2 r - 1, the h[(z - k)/2 r ] assume 
nonzero values on disjoint sets. On the other hand, h[(z — k)/2 T ] assigns 
different values to the points k and k + 2 r while any h[ (z — k )/2 r } assigns 
the same value to these points for r < r. Thus, h [(z — k)/2 r ] is not a lin- 
ear combination of other h[(z — k)/2"] fors ^ r. This establishes the in- 
dependence of the h's. Note that for ^ k < 2 T and ^ r < n, the func- 
tion h[(z — k)/2 r ] has period 2" and mean zero. By taking suitable 
linear combinations of the 2" — 1 independent h[(z — k)/2 T ], we can 
form functions / which assume any desired values on the points 0, 1, 2, 
• • • ,2" — 2. Of course, we must have 

2»— 2 

f(2 n - i) = - E/c*). 

z=0 

Consequently the 2" — 1 functions h[(z — k)/2 r ] form a basis for the set 
of all periodic functions with period 2" and mean zero. That is, any func- 
tion / with period 2" and mean zero can be written as a linear combina- 
tion of primitive functions with period 2". This completes the proof of 
(3). 
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To conclude this section we give an example which illustrates the ease 
with which the primitive components of a function may be found. Con- 
sider the function g given by: 



g: • • • , 1,2, — 5,7r,2,l — t,2,0, 



g has period 8 and mean zero. The only component h[(z — k)/2 r ] which 
can cause a difference in 0(0) and #(4) is 

fc(*/4): ••• ,1,0,0,0,-1,0,0,0, ••• . 

Since a h(z/±) assigns the points and 4 values which differ by 2a and 

g(0) - p(4) = \ 

then by choosing a = J we obtain the h(z/4:) component of g. Perform- 
ing similar calculations for h[(z — fc)/4], k = 1, 2, 3, we obtain 



4 V4/ V 2 / V 4 



+H^)-i*(^) 

given by 

(Ji (z) : 



33 — 7T 9 f 3 3 — T 9 7T 

>4>— 2 ' "4' 2'4'~ 2 ' ~4'2' 



(which has period 4). We apply the same arguments to the decomposi- 
tion of 0i (2) into the h[(z — k)/2], k = 0,1, and find 

given by 



M -? 5 -5 ? _ 3 3 _ 3 3 

0*W : ••' 4' 4' 4'4' 4' 4' 4' 4'"' 

so that g->(z) = —%h(z). Consequently g has been decomposed into 
primitive functions. Graphically we have: 
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We are ready to proceed to the n-dimensional generalizations of the 
results of Section II. The proofs given will use basically the same ideas 
as before although the technical details become somewhat more formal 
and involved. We begin with some definitions. 

Let Z n denote the ring of n-tuples of integers with component- wise 
addition and multiplication. That is, if a = (ai , ■•• ,a„) and b = 
(bi , ■ • • , b„) are elements of Z" then 



a + b = (fli + 6i , • ■ • , On + b n ) 



and 



a-b = (ai-bi , 



&„). 



In general, unless otherwise noted, lower case Latin letters without sub- 
scripts will denote elements of Z n ; lower case letters with subscripts will 
denote elements of Z, i.e., integers. If a e Z and q = (qi , qi , • • • , q n ) 
e Z n then we define aq to be {aq x , aq 2 , • • ■ , aq„). The w-tuple (1,1, 
• • ■ , 1 ) will be denoted by e. By a < b we mean a, < fr, for 1 ^ i ^ n. 
A function /: Z" — > R (the real numbers) is said to be 'primitive if 
there exist a,x (1) , ■■■ , x {n) e Z n such that 
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■_ l y 1 +-+c n j f z = CiX U) _|_ . . . + Cn . c <») + a 



otherwise 



for all z e Z n . 

Let g denote the real vector space generated by the set of all primitive 
functions. Z" will denote the subset of Z" consisting of those w-tuples 
which have all -positive coordinates. If m e Z n+ then P m is defined by 

P m = \xzZ n : ^ x < m) 

(i.e., ^ Xi < rrii for 1 f£ i ^ n, where will be used to designate both 
an element of Z and also the n-tuple (0, 0, • • • ,0).) 
A function /: Z" — > R is said to have -period m if 

f(z + km) = f(z) for all z, k e Z". 

If / has period m and 

scP m 

then / is said to have mean zero. Let 5 m denote the real vector space of 
all functions of period m which have mean zero. Next, we define 

f[(z ~ a)/b], b ^ 0, 

by 

' (0 otherwise 

For a e Z, a ^ 0, let E(a) and 0(a) denote the "even part" and "odd 
part" of a respectively. In other words, if a = 2 ,3 (2 i u + 1) for /3, n e Z 
then E(a) = 2* and 0(a) = 2/i + 1. For m = (m, , • • • , m„) e Z", 
E(m) will denote the n-tuple (E(mi), • • • , i?(ra„)) with 0(m) defined 
similarly. 

Finally, for m e Z", let 5 m * denote the real vector space generated by 
the set of functions 



\0(m)) : 



f [ow) :fe * Elm) ' aeZ 

We note that if m = e = (1,1, • • • , 1 ) then & m * = JF m ; in general, we 
always have fF OT * C fF„, . 

We come now to the main result of the paper. This is the following: 
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Theorem. 

8 (1 S m = £F m * 
for m e Z n+ . 

The proof of this theorem will proceed in a series of Lemmas parallel- 
ing the steps taken in the preceding section. 

Lemma 1. If g(z) eg then g[(z — a)/r] e g for all a e Z n and r e Z" + . 

Proof. We first show that if f(z) is primitive then/[(z — a)/r] is primi- 
tive. If we assume f(z) is primitive then by definition there exist 
b,x ll \ ••• , x M e Z n such that 

f (_!)«!+•••+«» if z = ClX w + . . . + Cn x {n) + b. 
[0 otherwise 
On the other hand 

(f(y) if z = ry + a 



(r?) 



r'~" 



f 

[0 otherwise 

Therefore 

( _ 1 )c 1+ ...+e B if z = r ( CiX w + . . . + CnX M + 6 ) + a 

otherwise 

f(_ i) fi i+- "+«» if 2 = Cl ( ra (1) ) + ■•.+c„(ra: (w) )+r6 + a 

[0 otherwise 

and hence, f[(z — a)/r] is primitive. By applying this result to a linear 
combination of primitive functions, i.e., an element of g, the lemma fol- 
lows. 

Let i,r: Z n — > R be defined by 

i(z) = 1, tt(z) = ziz 2 • • ■ z n for all z = (21 , z 2 , • • • , z») e Z". 

Lemma 2. Let c,m e Z n+ , r = cm and suppose f: Z" — » 7? /mis period 
m. Then for any a e Z n 
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Proof. We first note that since 

r = on = E(c)0(c)E(m)0(m) 

= E(c)E(m)0{c)0{m) = E(r)0(r) 
then 

E(r) = E(c)E(m) and 0(r) = 0(c)0(m). 
By definition 



s ' L L1 )h*> -s/w 

?of e) m \0(0(c)m)J £i 



tePoic) 



and 



where 



•Ml \0(W) / z B B 



A = {2: ^ z = W(c)0(m) + a < 0(c)m for some k\ 
and 

5 = {z: ^ z = kO(m) + a < m for some fc}. 

Hence, for each set, the values which k may assume are just a transla- 
tion of P«(m) , there being ir(2?(m)) values in all. Since 0(m) and 
0(c)0(m) are odd (i.e., each component is odd), then A and B both 
contain a complete residue system modulo E{m). Consequently, since all 
the elements of A and B are congruent to a modulo 0(m) then modulo 
E(m)0{m),A and B are identical. Since m = E{m)0{m) and / has period 
m then the sums S/( 2 ) an d ]C/( 2 ) are equal. 

zeA ztB 

We also note in general that for any s e Z n+ 

,2<wV w -* (£(c)) £<w)' fe) 

since Pb< C )„ is the disjoint union of w(E(c)) copies of P. . Therefore we 
have 



(5) 
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and the lemma is proved. 

We should note that as a corollary to this lemma we obtain : 

i ( J / (z) has mean zero iff 

\ cm / 

i ( ) f(z) has mean zero. 

\ m / 

We are now in a position to prove the important 

Lemma 3. Suppose g e g /ias period p e Z n+ . TVien for all a e Z", 

S'te)' w - ft 

Proof. We first show that the above conclusion holds if we assume that 
gr = /is primitive with period q = (a, a, • • • , a). In this case there exists 
c,x ll) , • • • , x {n) e Z n such that 

f (_l)«i+-+«n ^ 2 = 0ia .U> + .. . + fl||X <»> + c 

/CO = L . 

(0 otherwise 

To each u = (w x , • • • , u„) e Z n we can associate the unique point 
» = 0>i , • • • i v n ) e Z* suchthat^i = «i± #(?), w< = w» fori > 1, where 
the ± sign is chosen so that Z" is decomposed into the union of disjoint 
pairs { u,v} . It follows at once that 

Uyx {1) + • • • + u n x {n) m Vl x {1) + • • • + v n x {n) (mod 0(q)). 

Since 

t=i i=i 

is odd then 
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f( Ul x w + • • • + u n x (n) + c) - -ffaP + • • • + »„x (n) + c). 
Also, note that for all a e Z" 

l \ m ) = l \ m /■ 

Since / has period q then we must have 

as asserted. 

We may now remove the restriction that / has period of the form 
q = (a, a, • • • , a) = ae. If we assume / has an arbitrary period p e Z n+ 
then it is certainly true that / also has period ir(p)e. By above we have 

pf£u \0(T(p)e)/ 
Since p divides ir(p)e then by (5) we see that 

^te?)'^ - (6) 

Finally, to prove the lemma assume that 

t 

g = 2 <*jfi 

where the a> are real and the /,- are primitive. If /,• has period p U) then 
by (6) we have 

If # has period p and g denotes p (1) p <2) • ■ ■ p w p then by Lemma 2 

^tef)'' (z) = for 1SJ ' = ( - 

Therefore 

S<w)' w -S*S'(w)'' w -° 

so that applying Lemma 2 again we obtain 
since y has period p. This proves the lemma. 
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The final lemma is an n-dimensional generalization of (3). Its proof, 
however, is considerably more complicated. 

Lemma 4- If m = $Fe for some a e Z + then JF m C S- 

Proof. It will be sufficient to show that there exist 2" a — 1 functions 
in Q which also belong to 5 m and which are linearly independent over R. 
Let C lk) denote the k X k matrix of the form 



C (k) = 



1111 
1-111 
11-11 
111-1 



[1111- 

and let D {k) denote the k X k matrix of the form 



,(k) 



D w = 



In other words, 



and 



fl 



for i = 1 



Ca = 



da = 



1 - 25, 7 for i > 1 



f0 if i +j ^ k 
1 1 otherwise 



1 
1 
1 

1 

-1 


















1 














1 


1 











1 


1 


1 
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where 8 tj is the Kronecker S-function. Define B K W to be the v X v matrix 
of the form 

/o c {v - K+1) \ 

B K M = for 1 ^ k ^ v 

\D { - 1) / 

where denotes the appropriate zero matrix. Let r K ,\ " denote the point 
of Z" formed from the Xth row of B, w . Finally, let/,/* denote the func- 
tion in S defined by 

f(-l)" ,+ - + " n if z = a ir J v) + •• • + cr,., 00 

[0 otherwise 

We show first that the functions//"', 1 ^ k ^ v, are linearly independ- 
ent over R. To accomplish this it suffices to show that for any k, 1 ^ k ^ 
v, there are two points p and q in Z" such that 

while 

/r ( "(p) =/r W (e) for k<t£v. 

What we show in fact is that if |/, w (p) \ = 1 then f K+ i y) (p) = 1 for 
1 S* * < v. This may be proved by showing that if s e Z" is any Z-linear 
combination of the r K ,^ v) , 1 ^ X ^ v, then s can be written as a Z-linear 
combination of the r.+LX 1 "', 1 ^ X ^ p, such that the sum of the coeffi- 
cients is divisible by 2. We proceed by induction on v. For v — 2 we 
have 



Since 



and 



„ <2) _ „ (2) l v <2) 

ri.i = r 2 .i -+- r 2i2 



(2) (2) (2) 

n.2 = r 2 , 2 - r 2 ,i 



then any Z-linear combination of the ri,\ (2) can be written as a Z-linear 
combination of the r 2 ,x (2) with the sum of the coefficients divisible by 2 
and the assertion is true for this case. Now assume the hypothesis for v 
and let 



s = £ axr,, 



(-+D 
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E 

be a Z-linear combination of the r K ,\ v+l) . There are two cases. 

(i) Suppse k ^ 2. Note that B K l ' +1) and B K+1 l ' +1) have the following 
forms: 



B K 



(h-u 




11 ••• 10 •■• 



K - 1 



(.•+1) _ 







We also note the important fact that r lc+ i,„_ )[ ( '' +1) has all zero components 
except for the xth component which is 1 . Thus we have 

'iC.V+l 'K+l.f+1 ' H+l.K— K 

Now s can be written as 



Ed'+D I 



x=i 



+1 



("+D 



By the induction hypothesis (since the r Ki \' ,+l) differ from the r^x*" 5 , 
1 ^ X ^ j/, only in an extra zero component) we can write 



ZW +1) = X>xr. +1 .x ( 

x-i \=1 



y+l) 
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where ^ b\ is an even integer. Hence 
x=i 

V> , (»+n , „ / (»-+i) „ (>+\)\ 



x=i 



Here, s has been written as a Z-linear combination of the r« + j,x " with 
the sum of the coefficients even and consequently this case is finished, 
(ii) Suppose k - 1. We write Bi" +n and B 2 ( ' +1> as 



Bi 



(••+1) 






1 


1 


... i 


1 - 


1 












1 












1 


B 2 ( " +1) = 






B 2 M 




1 








1 


... i 


1 - 


-1 



We note that 

As before, by the induction hypothesis we can write 

L(»+D Y^ L „. ("+D 

axn.x = 2^ Oxr 2 .x 



where £] &x is even. Consequently we have 
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x=i 

— V h v ( " +1) _l_ r, (m ( "+ 1) _L *■ ( ' +1) \ 

X-l 

which expresses s as a Z-linear combination of the r 2 ,x (,,+l) with an even 
coefficient sum. This completes the induction step and the proof of the 
assertion that the/« , 1 ^ k ^ ?, are linearly independent over R . 
A more careful examination of B ie ( ' ) reveals the following: 

(a) r K ,x + i W - r«.x w = (5 liV _x , 8 2 ,_x , • • • , 8,,_x) 

(b) r«,i ' — r«,x " = (25i,, + x_i , 25 2>lt+ x_i , • ■ • , 25„, (t+ x-i) 
for 2 ^ X ^ k - A : + 1. 

(c) Er,/" - (v-k - 2)r, a w = (25 li)t , 28, .« , ••• , 2*,.,). 

X=2 

Since the linear combinations of the r K ,\ r) in (a), (b) and (c) all have 
the sum of coefficients an even integer then 

/. w «Mu,%,x,--' ,2fi,, x )) =f. iv) (0) 

for 1 ^ X, k ^ v. Hence /« has period 2e = (2, 2, • • • , 2). Also we note 
that the only points in F 2c (Z Z" at which /, is nonzero are just those 
Z-linear combinations of the r«,x which have all coordinates or 1. It 
is not difficult to see that the only points of this type which may be gener- 
ated are the 2* points of the form (ci , c 2 , • ■ • , c,_i , c , c , • • • , c ) 
where c, = or 1. By a translation of f K (v) by a we mean the function 
fj v) denned by 

fj v \z) = /,«(• - a). 

By letting the a range over the set of points 

A K = 1(0,0, ••• ,0, rfi ,d 2 , • • • , d_»):<*x = or 1}, 

« 

the 2"~* translations f K , a " , a e A K , have the property that for each p e 
P 2e , exactly one of the/«, assumes a nonzero value at p. In fact, if we 
define an inner product (f«,o W , /x,6 W ) for/«,« w and/x.b 00 by 

pcP 2e 

then the inner product of any two distinct functions /«,o ,/x,& j a e i,, 
& e ^.x , 1 ^ k, X ^ j/, is zero. Since 

(fj*,fj*) = 2< 
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then by introducing the "normalized" basis functions 

Jk.o — — jK.a 

we see that the set of (1 + 2 + V 2"- 1 ) = 2' - 1 functions /,,«<->, 

ae A, , \ ^ x g v, are orthonormal. Thus, we have shown that F 2e a g. 
The extension of this technique to show that F Vr C g is quite similar 
and will be omitted. The basic idea is simply to introduce the "expan- 
sions"/,^ 00 <rf/«.. W defined by 

/«.«.x W (z) = JJ V) fofc) for X = 2°, 2\ • • • , 2"" 1 , 

and then by taking suitable normalized translations of these functions, 
obtain an orthonormal basis (in g) for SFa-e . This shows that 

s„, = 5 2 a c c g 

and the proof of the lemma is completed. 

We are now ready to proceed to the proof of the 

Theorem. 

g fl SF m = QF»* 

for m e Z" + . 

Proof: g H J,,, cz $,„ * 

Let / e g D g r „, . Since / has period m then by Lemma 3, we have for all 
a e Z" 



£ \55S?) 



/(*) =o. 

Since 



' w -JS.,te> w 



and each of the functions i[(« - a)/0(m)]f(z) can be written as 
h[(z — a)/0(m)] for some h e 5 B ( m ) then/ e 5 m * and this direction is 
established. 

griff,,, =3 ff,„* 



We have already noted that "S m * C ff m . It remains to show that ?„,* CI g. 
By definition SF»* is the real vector space generated by the set 

j/i[(2 - a)/0(m)]: h e ffs (m) , a e Z n } . 
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By Lemma 4 we have 

Thus, if h e ??*(„,) then /i e 9. But by Lemma 1,/ieg implies 

fc[(2-a)/0(m)]e8. 

Therefore, since 9 contains a set of generators for ?„,* then 5 m * CZ 9. 
This completes the proof of the theorem. 

IV. CONCLUDING REMARKS 

As a concluding example of the results of the preceding section, we 
consider the decomposition of the function / generated by the charge dis- 
tribution of the crystal structure of potassium tantalate, KTaOs . This 
compound forms face-centered cubic crystals with a charge distribution 
as shown in Fig. 1. That is, a +1 is situated at each vertex, a —2 at each 
face-center and a +5 is located in the center of the cube. The periodic 
function / defined by this distribution has period (2,2,2) and is shown in 
Fig. 2 (which is the forward upper left octant of Fig. 1 ). We have 



Bx (3) = 



1 


1 


1 


) 


1 


1 


) 


1 


1 


-1 


1 


, B^ = 


1 


-1 


R (3) 


1 


1 


1 


-1 




1 







1 1 



so that the 7 basis functions into which / will be decomposed are as shown 
in Fig. 3. 




Fig. 1 — Charge distribution of KTa0 3 . 



LATTICE-PERIODIC FUNCTIONS 



1213 




Fig. 2 — A period of the periodic function. 

-I -I 



/ 














y 



1 i, (o,o, o) 



+10 



-I 



/ 






( 



f (3) , ^ 

l, (1,0,0) 



+ 1 

-I 



I, (0,1,0) 



f (3l 
' 1, Utii< 



+1 

-if*- 



f l3 \ 1 
1 2,(0,0,0) 



-1 



1 2,(0,1,0) ' 3,(0,0,0) 



Fig. 3 — The seven basis functions. 
+ 1-1=0 -1-1 =-2 



-2 + 2 + 1 = 1 



-1-2 + 1 =-2 — - 



-+1-2 + 1 = 



+2+2+1 = +5 



+ 1-1=0 -1-1=-2 

Fig. 4 — Decomposition of the periodic function. 
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The coefficient of /i,<o,o,0) is obviously \ (0 — ( — 2)) = 1, etc., so 
that we obtain 

t r (3) r (3) r (3) 

/ — Jl, (0,0,0) /l, (1,0,0) /l, (0,1,0) 

T 4/1.(1.1.0) — 4/2.(0,1,0) — /3, (0,0,0) 

Graphically, this equality is shown in Fig. 4. 

The author gratefully acknowledges many enlightening discussions on 
this subject with H. O. Pollak (whose ideas along the lines of generating 
functions led to a short solution of the one-dimensional problem) and 
W. J. C. Grant (who originated the problem). 
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